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Abstract 

We compute the leading QCD corrections to K-K mixing in the supersym- 
metric standard model with general soft supersymmetry-breaking parameters. 
We construct the /S.S = 2 effective Lagrangian for three hierarchies of super- 
symmetric particle masses, namely, when the gluino mass is comparable to, 
much greater than, or much less than the masses of the first two generation 
squarks. We find that the QCD corrections tighten the limits on squark mass 
splittings by more than a factor of two. 
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1 Introduction 



Low-energy supersymmetry is a leading candidate for physics beyond the standard 
model because it stabilizes the gauge hierarchy. However, the minimal supersym- 
metric extension to the standard model contains over 100 new parameters. Present 
experimental searches are beginning to place significant constraints on the parameter 
space. To date, the most important limits arise from the study of fiavor-changing 
neutral currents (FCNC). 

In general, supersymmetric particles give rise to large contributions to fiavor- 
changing neutral currents Jll. The present limits on these processes place strict 
constraints on the squark and slepton mass matrices. Indeed, the discovery of su- 
persymmetric FCNC will be an important step towards understanding the source of 
supersymmetry breaking. 

There are several FCNC processes which constrain the parameters of supersym- 
metric models. For example, the rare decay 6 — 57 has been used to constrain the 
mass parameters of the third generation squarks [0]. In this paper, we will focus on 
K-K mixing, which gives the most stringent limit on squark masses of the first two 
generations. 

In what follows, we will extend the analysis of Refs. to include the lead- 

ing order QCD corrections to K-K mixing in the minimal supersymmetric standard 
model (MSSM). The QCD corrections are clearly very important, and for the stan- 
dard model, complete leading order and next-to-leading order analyses have 
been carried out. (Our results can be easily generahzed to B-B and D-D mixing, in 
which case they also involve the third generation squarks.) 

For general squark masses, the most important supersymmetric contribution to 
K-K mixing comes from the gluino box diagrams.| We will, therefore, restrict our 
attention to the gluino diagrams, and ignore all diagrams with chargino and neu- 
tralino exchanges. These diagrams depend on the details of the Higgs sector and are 
generally less important than the diagrams we consider here. 

As is by now a standard procedure, we define the source of FCNC in terms of 
dimensionless fiavor-changing insertions 5, which parametrize small deviations from 
the case of fiavor-diagonal soft squark masses. We follow the notation of Ref. but 
we omit the generation and weak isospin labels because we concentrate solely on the 
first two generations. 

This paper is organized as follows. In Sec. |^ we review the general formalism for 
calculating QCD corrections in the context of effective field theory. We then consider 
three cases for the hierarchy of the superparticle masses: (a) Mg ~ Mg, where the 
gluino mass, M^, is on the order of the average squark mass, M^, as well as (b) 
Mg -C Mg and (c) Mg ^ Mq. In each case, we construct the effective Lagrangian at 
the matching scale and compute the renormalization group equations for the Wilson 

■'•In the minimal supergravity model, the chargino and neutralino box diagrams can be of com- 
parable importance . For a recent review, see Ref. . 
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coefficients. In Sec. |^ we present the results of our numerical analysis. We find that 
the QCD corrections tighten the limits on the squark mass splittings by more than a 
factor of two. (Of course, our results are subject to the usual uncertainties associated 
with the hadronic matrix elements.) We reserve Sec. § for our conclusions. 



2 QCD corrections to the effective Lagrangian 

2.1 General formalism and conventions 

To analyze a low energy physical process in a theory with several mass scales, it is use- 
ful to construct an effective field theory, in which particles with masses greater than 
the scale of interest are integrated out . This gives rise to an effective Lagrangian, 
which can be written as follows, 

/:eff = ECA(/i)OA(/i), (1) 
A 

where /i is the renormalization scale, and the operators Oa involve only low energy 
fields. The Wilson coefficients Ca are obtained by matching S-matrix elements in 
the full and effective theories at the threshold (matching) scale. The heavy particles 
modify the ultraviolet behavior of the low energy theory, but the infrared physics is 
the same in both. 

In any realistic calculation, the Wilson coefficients and operators must be evolved 
from the matching scale to the scale associated with the low energy process under 
consideration. The renormalization group equations are given by0 

^dOA ^ ^ {^j't^^^^Cb, (2) 

dfi djj, 

where 7 is the anomalous dimension matrix. Using the one-loop RGE for the strong 
coupling, 

(3) 

where 

di = + -n„ + -ns + 2n-A (4) 

167r2 V 3 ^ 6 * V ^ ' 

and the number of active quark, squark and gluino flavors, respectively, 

one can find a formal solution for the Wilson coefficients at any scale /z. 



C(/i) = exp 



C{M). (5) 



^We do our calculation in the MS-scheme. 
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Here M is the matching scale and 7 is the one-loop anomalous dimension matrix, 
which is of order g1. This procedure resums the leading logarithmic QCD corrections. 

In the following subsections, we will use the effective field theory formalism to 
construct the leading-logarithmic QCD-corrected — 2 effective Lagrangian in the 
low energy limit of the MSSM. We will consider models with the following hierarchies 
between the gluino and first two generation squark masses: Mg ~ M^, Mg -C Mq and 

M-g > Mg. 



2.2 M~g ~ Mq 

When Mg ^ M^, we choose to integrate out the squarks and gluino at the scale Msusy, 
which we define to be the geometric mean of the squark and gluino masses, AfsusY = 



MgMq. The supersymmetric contribution to the AS — 2 effective Lagrangian is 
then 



'^.«=^^'^^i:cmoaw, (6) 



where we have defined the operators 







02 






— 1^ 




— It J y J 


0, 





(7) 



plus other operators 0i,2,3, with the obvious exchanges L ^ R m. Ci,2,3- In these 
expressions, the superscripts i,j are SU{3) color indices. All other operators with 
the correct Lorentz and color structure can be related to these by operator identities 
and Fiertz rearrangements. 

At the matching scale MgusY, we determine the Wilson coefficients to be 

Ci(Msusy) = (24x/6(x) + 66/6(x)) 5l^, 

C2(MsusY) = 204a;/6(a;)4i, 

C3(MsusY) = -36x/6(x)4i, (8) 
C4(Msusy) = (504x/6(x) - 72/6(x)) SllSrr - lS2Mx) SlrSrl, 
C5(Msusy) = {2ixfe{x) + 120fe{x)) SllSrr - ISOMx) SlrSrl, 
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where x = M~/M? and the 5's come from insertions of the squark mass matrix. The 
functions /e and /e 

, , , 6(1 + 3a;) log x + ~ 9x'^ - 9x + 17 
fe{x) = 



6(x- 1; 
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~, , 6x(l + x) log X — s'^ — 9a;2 + 9a; + 1 

= 3(^3Tji (9) 

arise from momentum integrals in the gluino-squark box diagrams. They have the 
limits 

a;/6(a;) ^ 0{x), U{x) ^ (x « 1) 

:r/6(x)^-^, Mx)^Oix~^) (x»l), (10) 
ox 

so the leading behavior of the effective Lagrangian goes as 1/M| for x <^ 1 and 1/M| 
for X ^ 1. The values of the coefficients Ci 2,3 are obtained by replacing L ^ R 
in eq. (^. Note that our result disagrees with that of Ref. 0], but confirms that of 
Ref. §. 



The operator Ci is of — A type. Its anomalous dimension is well known [|C 

liO,) = (11) 
The one-loop anomalous dimensions for the other operators are 

^(O^Os) = ]], (12) 



9? 


-r: 


127r2 


4 8 


9! ( 


-8 


87r2 I 


-3 1 



^ ( -3 1 j ■ 

With these anomalous dimensions, we evolve the Wilson coefficients to the scale //had, 
where the hadronic observables are defined. We find 



Cl(/^had) 


= 77iCi(Msusy) 




C2(/^had) 


= 7722^2 (^SUSY 


) +'723C3(^SUSy), 




= ''?32C2(^SUSY 


) +^33C3(MsuSy), 






+ ^{Vi - V5)C5{MsvSY 









(14) 
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where 



^ \a.,(/iha,d)/ \a,,(m.)/ \aJmf,)J \ aJn 



6/21 

— ) 

^(/ihad)/ \as{mc)J \as{mb)J V as("^t) / 
ri22 = 0.983r/2 + O.OlTr^a, V23 = -0.258ri2 + 0.258ri3, 
r/32 = -0.0QAr]2 + 0.064r/3, r/33 = O.OlTr^s + 0.983r/3, 

-2.42 2.75 -4 1/2 /Tr-N 

V2=Vl , ^3 = ^1 , ^4 = , V5=Vl ■ (15) 



2.3 <C (x <C 1) 

When Mg <C Mg (x <^ 1), we proceed in two steps. We first integrate out the heavy 
squarks at their own mass scale, Mg (see Fig. |l|). This gives an effective Lagrangian 
with AS* = 1 and AS = 2 four-fermion operators. The AS = 1 operators are of the 
form 

^{T^T^^ng^fsl (16) 

together with similar operators with the obvious exchanges L ^ R and/or d s. 
The AS = 1 operators give rise to one-loop AS* = 2 diagrams which are suppressed 
by M^/M^, so we ignore them in what follows. The AS* = 2 terms in the effective 
Lagrangian are given by 

Ces = ||^|-225i^Oi-224^ai + Wij«(24O4-40a) 

+5L/?5/?L(44O4 + 60a)}, (IT) 

where the matching is performed at the scale Mg. 

To construct the effective Lagrangian, these AS* = 2 operators must first be 
evolved to the gluino scale, and then to the hadronic scale.' The one-loop anomalous 
dimensions for operators ^1,4,5 and Oi were calculated previously, so the leading 
AS" = 2 effective Lagrangian at the hadronic scale can be obtained by scaling eq. (|T^). 
We find 

= "^'(^g) |_22^2 ^ _ 225lr,K, Oi 
216M|l LL 1 1 RR 1 1 

g 

+SllSrR (-(4k4 + 5/€5)C4 - 4OK5C5) 

^Integrating out the gluinos does not induce any new AS — 2 operators. However, it changes 
the beta function for the strong coupHng. We have assumed that the third generation squarks have 
masses of order Mg, in which case the one-loop beta function coefficient is —13/3 between Mg and 
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AS = 1 
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O 



Ml 



+ 



A5 = 2 ^ ^ 
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O 



Ml 
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Figure 1: The matching of the S-matrices in the fuh and effective theories, at the 
scale Mq, for the case of Mg -C Mq. At one loop, the AS* = 1 operator generates a 
subleading AS* = 2 operator, of order M?/M^. 



where 



Ois (/^had 
4 



+Slr6rl ((64«:4 - 20^5)04 + 60/C5O5)|, (18) 

6/27 / (m,) ^^ 6/25 ^ (^^) ^ 6/23 ^ a^(Mg) ^ 6/21 ^ ( jy^,) ^ 6/13 

(19) 



asirric)/ \as{mb)J \as{mt) J \as{Mg) 



_ 1/2 
^5 — '^1 • 



2.4 > Mq {x > 1) 

When Mg ^ Mg (x ^ 1), we must again proceed in two steps. We first integrate 
out the gluino at Mg (see Fig. |^). The effective Lagrangian between Mg and Mg, 

C,s = T.DAif^)QA{l^), (20) 

A 

contains the following leading order operators 

Qs = n^csr, (21) 
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AS = 




AS = 2 



I I 
I I 




o 



o 



Figure 2: The matching of the S-matrices in the fuh and effective theories, at the 
scale fi = Mg, for the case Mg S> Mg. The AS" = operator generates a leading 
AS = 2 operator, whose coefficient exactly matches the a; S> 1 limit of the corre- 
sponding box diagram. This implies that the Wilson coefficient of the first AS' = 2 
operator on the right-hand side is zero. 



as well as other operators with obvious exchanges, L R and/or s <-> d. The 
superscript c in Q3 stands for the charge conjugated field. 

By matching at the scale Mg, we determine the Wilson coefficients to be 

D, = -3L>2 = GDs = (22) 

Mg 

The one-loop anomalous dimensions are 

o2 / -8 \ 3g2 
7(QiQ2) = I, _J, 7(Q3) = -^. (23) 

Therefore the Wilson coefficients Da (A = 1, 2, 3) at the scale Mg are given by: 

D,(Mg) = SiD,(M~g) + ^(ei-e2)D2(M~g) 

DA{Mg) = SADA{Mg), A = 2, 3 (24) 

where 



= ' = ' = • ^''^ 
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We now evolve the effective Lagrangian to the scale Mq, where we integrate out 
all the squarks. At that scale, the AS" = 2 effective Lagrangian is given by 



eff 



216M| 



+ 



(26) 



This effective Lagrangian matches the Mq <^ Mg limit of the full theory. 
Finally, this Lagrangian must be evolved to the hadronic scale. We find 



eff 



216M| 



+ 



4l( (^(64e? - el) - 8el)r^',, + {2sl - 8el)r^',,)02 



(3(645^ - e',) - 8ei)v',2 + i^e^ - 8ei)v',, )Os 



+Ae',r][6iL0, + {L^R,0^0) 



(27) 



where the ri'^s have the same form as the ?7's, with the scale Msusy replaced by Mq 
in eq. (plSf) . 



3 Numerical results 

In this section, we will use the AS" = 2 effective Lagrangians to compute the super- 
symmetric corrections to the i^-meson mass difference, 

ArriK = 2Re{K\Ces\K). (28) 

We shall see that the measured value, Atjik = 3.5 x 10^^^ MeV, places tight limits 
on the parameters 6. 

Note that the AS* = 2 effective Lagrangians can also be used to compute the 
supersymmetric contribution to the CP- violating parameter ex, 

^ lMK\C..m .29) 

If the phases of the (5's are of order one, the experimental limits on ex place very 
tight constraints on the squark mass matrices (see, for example, Ref. [||). 
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X 




Vl Re 52^1 


0.3 


2.3 X 10-2 


3.0 X 10-2 


5.5 X 10-3 


2.8 X 10-3 


1.0 


5.0 X 10-2 


6.6 X 10-2 


6.3 X 10-3 


3.2 X 10-3 


4.0 


0.12 


0.16 


9.2 X 10-3 


4.6 X 10-3 












0.3 


3.1 X 10-^ 


1.0 X 10-3 


4.2 X 10-3 


1.4 X 10-3 


1.0 


3.6 X 10-3 


1.2 X 10-3 


7.2 X 10-3 


2.4 X 10-3 


4.0 


5.3 X 10-3 


1.7 X 10-3 


1.7 X 10-2 


5.6 X 10-3 



Table 1: Limits on the Re^'s, for squark mass Mg = 500 GeV and different x = 
M?/M~. In each case, the left (right) numbers are limits without (with) leading 
order QCD corrections. 



In each of these expressions, the effective Lagrangian jCeS contains the operators 
Oi through O5, with coefficients evahiatcd at the hadronic scale, which we define to 
be the scale where a, = 1. We evaluate the hadronic matrix elements and find 



{K\0,\K) 
{K\0,\K) 
{K\Os\K) 
{K\04K) 
{K\0,\K) 



^rriKfKBi, 



5 



24 \ms + rud 
1 / rriK ^ 2 



24 Knis + rrid 



1 

24 



rriK 



rus + md 
rriK 



1 1 

8 12 KrUg + md 



mKfKB2, 



(30) 



The coefficients Si_5 characterize the long-distance hadronic physics; we take Bi = 1, 
as determined in the vacuum insertion approximation. 

In the rest of this section, we use these results to illustrate the size of the leading 
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\ \ I \ \ I I \ I \ 

1 10 100 

M~ [TeV] 

Figure 3: Limits on the Re(5's as a function of the common first two generation 
squark mass, Mq, for a hght gluino mass of Mg = 200 GeV. The sohd fines describe 
our effective field theory result. The dashed lines correspond to decoupling the 
supersymmetric particles at MsusY = i/M^Mg, without including the leading order 
QCD corrections. 



order QCD corrections to K-K mixing, for each of the three cases discussed above. 
For our numerical work, we take fx = 160 MeV, rrix = 498 MeV |TT| and the current 
mass rUs = 150 MeV.|||] We use the strong coupling as determined from electroweak 



measurements, as{Mz) = 0.118 |jTT[. This gives ^^(mc) = 0.35 and ^^(mf,) = 0.22. 

We first consider the case Mg ~ Mg. The QCD-corrected AS" = 2 effective 
Lagrangian gives rise to the limits shown in Table 1. In the table, we have taken 
the first two generation squark mass Mq = 500 GeV and varied x = M^/M~. The 
left (right) numbers correspond to the limits without (with) the leading order QCD 
corrections.^ 



The most stringent limit in the table comes from y | Ile6LL6jiji\. For a given mass 
spUtting, we see the QCD corrections increase the squark lower bound by a factor 
of three! Note that the corrections are generally more than 70%, so they cannot be 
ignored in computing limits on the squark masses. 

We next turn to models with very different squark and gluino masses. In Fig. |^, we 

"This mass is not well determined; see Ref. jl^ for recent analyses. Decreasing nis to 100 MeV 
tightens the bounds by at most 30%. 

**Our numbers agree with those of Ref. |^ if we foUow their procedure and decouple the squarks 
and gluino at Mz and neglect the QCD corrections. 
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T 




^ I \ I \ \ \ 

1 10 100 

[TeV] 

Figure 4: The same as Fig. |3|, as a function of Mg, for light first two generation 
squarks, Mq = 200 GeV. 

consider the case of a Ught gluino, and plot the hmits on 5's versus the squark mass. 
The sohd hnes describe our results, while the dashed lines correspond to decoupling 
the supersymmetric particles at Msusy = ^MgMg, with no QCD corrections. 

From the figure we see that the QCD corrections tighten the most stringent bounds 
by about 50%. For Mg = 200 GeV, with 6 of order 1, the bound on the squark mass 
is Mq > 200 TeV. Such a heavy squark can drive the third generation squark/slepton 
mass squared negative if supersymmetry breaking is transmitted by gravitational 
interactions |TB| . 

Alternatively, one can make the gluino heavy while keeping the squarks light. We 
show limits for this type of model in Fig. ^ For Mg = 200 GeV, and 6 of order 1, 
the bound on the gluino is Mg > 200 TeV. (Note, however, that in this scenario, the 
charginos and neutralinos must also be heavy to suppress FCNC's from electroweak 
diagrams.) In a unified model, a large gluino mass drives up the squark masses 
through their renormalization group evolution, so this scenario can only be natural 
in models with TeV-scale supersymmetry breaking. 

4 Conclusion 

In this paper we have calculated the leading order QCD corrections to supersymmetric 
contributions to K-K mixing. We find the corrections to be significant. Indeed, for 
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the case Mq ~ Mg, the QCD corrections increase the lower bound on the first two 
generation squark masses by a factor of three. 

For the case Mq » Mg, we find that QCD corrections also increase the squark 
lower bound. This exacerbates the naturalness problems associated with such a 
hierarchy, and can destabilize the effective potential when supersymmetry breaking 
is transmitted by supergravity interactions. Similar results hold when Mq <^ Mg. In 
each case, the QCD corrections are important and must be included when deriving 
constraints from FCNC processes. 

We would like to thank A. Falk and M. Booth for very helpful discussions, and 
K. Agashe and M. Gracsscr for finding a mistake in a previous version of the paper. 
This work was supported by the U.S. National Science Foundation, grant NSF-PHY- 
9404057. 
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